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Partial-wave representations of laser o 

beams for use in light-scattering calculations 

Gerard Gouesbet, James A. Lock, and Gerard Grehan 


In the framework of generalized Lorenz-Mie theory, laser beams are described by sets of beam-shape 
coefficients. The modified localized approximation to evaluate these coefficients for a focused Gaussian 
beam is presented. A new description of Gaussian beams, called standard beams, is introduced. A 
comparison is made between the values of the beam-shape coefficients in the framework of the localized 
approximation and the beam-shape coefficients of standard beams. This comparison leads to new 
insights concerning the electromagnetic description of laser beams. The relevance of our discussion is 
enhanced by a demonstration that the localized approximation provides a very satisfactory description of 
top-hat beams as well. 


1 . Introduction 

Many optical particle-sizing techniques rely on the 
interaction between laser beams and the particles 
that are being studied. Examples are provided by 
Gaussian laser beams used in phase-DoppIer instru- 
ments, 1-3 laser sheets used in particle-image velocim- 
etry, 4 - 6 and top-hat beams used in the so-called 
top-hat technique. r - l ° If the diameter of the spheri- 
cal particles that are being studied is comparable to 
the characteristic beam width, the theoretical analy- 
sis of the light-scattering signature of the particles 
must rely on generalized Lorenz-Mie theory | GLMT) 
rather than on the usual plane-wave Lorenz-Mie 
theory. 11 A background in GLMT and its applica- 
tions may be gained from Refs. 12-15. 

In this formalism, shaped beams such as laser 
beams are mathematically expanded in terms of 
partial waves. The complex number that describes 
the amplitude and the phase of each partial wave in 
the expansion is called a beam-shape coefficient (BSC). 
These coefficients may be expressed as angular inte- 
grals of the radial component of the beam’s electric 
and magnetic fields. 13 Unfortunately, none of the 
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commonly used mathematical descriptions of laser 
beams is an exact solution of Maxwell’s equations. l& ~ 18 
As a result, the BSC’s produced by the angular in- 
tegration of these fields retain a weak dependence on 
the radial coordinate, 19 in contrast with the fact that 
the derivation of the partial- wave decomposition de- 
mands that these coefficients be constants. The 
residual radial dependence is an artifact that results 
from the imperfect description of the electromagnetic 
fields of the beam. 

On the other hand, a surprisingly accurate approxi- 
mation to the BSC’s is the so-called localized approxi- 
mation. 20 It is a simple analytical expression whose 
accuracy is typically approximately 1 part in 10 5 from 
the value of the constant portion of the BSC’s ob- 
tained by numerical integration. 21 This approxima- 
tion relies on the localized interpretation of partial- 
wave expansions, 22 which is an analogy to van de 
Hulst’s localization principle in Lorenz-Mie theory. 23 
The localization approximation BSC’s axe constants 
as is required of partial- wave expansions, and the 
beam descriptions generated by insertion of these 
BSC’s into the partial- wave expansion have been 
termed localized beams. 19 Such beams exactly sat- 
isfy Maxwell’s equations because they are built from 
basis functions with constant coefficients. 

Until recently, there had been no rigorous justifica- 
tion of the validity of the localized approximation. 
Initially, its validity was demonstrated by a compari- 
son of the numerical values of the BSC’s evaluated by 
the localized approximation with the values obtained 
by numerical integration or the finite senes 
method. 2124 " 26 When a beam propagates along the z 
axis of the coordinate system used to describe the 
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partial-wave expansion, i.e., an on-axis beam, one of 
us recently succeeded in giving a derivation of the 
localized approximation 21 that relied on the station- 
ary phase method, in analogy with van de Hulst’s 
derivation of the localization principle. 23 The deriva- 
tion, however, could not be generalized to off-axis 
beams, i.e., beams that are propagatingparallel to but 
not along the 2 axis. By the use of another technique 
that relies on Taylor series expansions, a final deriva- 
tion of the validity of the localized approximation for 
a focused Gaussian beam has been recently obtained 
for both the on-axis 19 and the off-axis cases. 27 This 
derivation uncovered a modification of the localized 
approximation that has been called the modified 
localized approximation. 

A significant ingredient in the derivation was the 
discovery of the so-called standard-beam expressions. 
We strongly believe that standard beams will prove to 
be the best mathematical description of Gaussian 
beams. Standard beams exactly satisfy Maxwell's 
equations, because they are constructed from BSC’s 
that are also constants. For the on-axis case, stan- 
dard-beam BSC’s are given by a simple infinite se- 
ries 19 that results from an extrapolation of the Davis 
procedure 16 for description of the electromagnetic 
fields of a focused Gaussian beam. The standard- 
beam BSC’s for the off-axis case have not yet been 
discovered, to our knowledge. 

In this paper we consider two aspects of the partial- 
wave representation of laser beams for use in GLMT 
scattering calculations, (a) In the context of an 
on-axis focused Gaussian beam, we examine the 
convergence properties of the infinite series that 
describes the standard-beam BSC’s. We also com- 
pare the values of the localized approximation and the 
modified localized approximation analytical expres- 
sions for the BSC’s with the standard-beam BSC’s 
that we use as a benchmark. We claim that the 
closer the localized approximation BSC’s come to the 
standard-beam BSC’s, the more accurate the local- 
ized approximation is in describing a focused Gauss- 
ian beam, (b) We apply the localized interpretation 
of partial-wave expansions to a top-hat beam and 
assess the accuracy of the resulting top-hat-beam 
localized approximation. To avoid burying the es- 
sence of the physics in complicated mathematical 
expressions, only the on-axis case is considered in this 
paper. 

This paper is organized as follows. Section 2 
summarizes the Davis formulation for description of 
the Gaussian beam electromagnetic fields and intro- 
duces the localized and the standard beams. Section 
3 compares the numerical values of localized and 
standard-beam BSC’s, providing new insights as to 
the nature of Gaussian beams and leading to the 
conclusion that standard beams should indeed be 
taken as the very definition of Gaussian beams. 
Section 4 further supports the validity of the localized 
interpretation of partial-wave analyses by describing 
the building of localized beams that provide a very 
satisfactory description of top-hat beams. 


2. Davis, Standard, and Localized Descriptions of 
On-Axis Focused Gaussian Beams 

A. Davis Formulation 

A description of the electromagnetic fields of a fo- 
cused Gaussian laser beam is provided by the Davis 
formulation. 16 * 1T<28 We consider a Gaussian beam 
that is propagating along the z axis from negative z' 
to positive z' (Fig. 1). Two parallel Cartesian coordi- 
nate systems must be used in this problem: (a) x', 

y', z', which is attached to the Gaussian beam and 
whose origin is at the center of the beam waist, and 
(b) x, y, z, which is used to describe the partial- wave 
expansion of the Gaussian beam. The origin of the 
x'y'z’ system with respect to the xyz system is z 0 . 
We start by considering the simplest case, z 0 = 0, and 
below we examine the more general on-axis case with 
z 0 * 0 when appropriate. 

We consider a monochromatic light wave with an 
exp(4*tio/) time dependence. This time dependence 
will be omitted hereafter, as is the normal practice. 
In the Davis formulation a laser beam is described by 
a linearly polarized vector potential, 


A = (A I , 0, 0). 

(i) 

The nonzero component A, is given by 


A x = \|/(x, y, z)exp( -ikz). 

(2) 





Fig. 1. Two coordinate systems that describe a focused Gaussian 
beam that is propagating along the z axis. The origin of the xyz 
coordinate system is at the center of the beam waist, and the 
partial-wave expansion is carried out with respect to the x, y , z 
coordinate system. 
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The function tMx t y , z) is unknown and must be 
determined. Such a determination will involve spa- 
tial derivatives. However, the transverse coordi- 
nates x and y scale with a small transverse character- 
istic length iu 0 , and the coordinate z scales with a 
large longitudinal characteristic length /. The scal- 
ing lengths Wq and l are taken to be the beam- waist 
radius and the spreading (or diffraction) length kw 0 2 t 
respectively. Rescaled dimensionless coordinates 
(€, t|, 0 may therefore be introduced according to 


Eq. (10), the fundamental mode ib 0 alone completely 
determines the vector potential A from which electric 
and magnetic fields are derived by the use of 




E 

= — V(V • A) - 

iu>A, 

fin 



H 

= (V x A)/>, 


(12) 

leading to 






■ 

J 

3 2 <Jio \ 



H 

ji 0 


llM + 7 + • ■ • 

-del \ 

exp(-ifcz), 

(13) 


The rescaled spatial derivatives Mi/dri, and 

M/ are now of the same magnitude. 

Within the Lorentz gauge, the vector potential A 
must satisfy the Helmholtz equation, 

V 2 A + k 2 A = 0, (4) 





d 2 ilJo 


+ s 4 


a 2 ii/ 2 

a^ 


. exp(-i^ 2 ), (14) 


E z ~ Eo 


-is 


dillp 


Jdx Jj 2 . a 2 ^ 0 '\ 

+ 1 


expi -ikz). 


providing the partial differential equation for iji: 


a 2 








5 + Hi* ~ + 777 = 0 






(5) 


In Eq. (5) we have introduced the small dimensionless 
parameter s given by 


S = a/g/l = l/kw 0 . 


( 6 ) 


H x = 0 , 


H y = 

H 5 = H 0 


, ,| , . *l»o 

U»o + S-| ll»2 -f £ — 
»(> 


. 5lll 0 
-is — 

dT] 


- is 3 


3(J/o 

df| 


(15) 

116) 

. . exp( -ikz), 1 17) 

]exp( -ikz). (18) 


Because $ is the ratio of two characteristic length 
scales that define the overall aspect of the beam, we 
name it the beam-confinement factor. For a plane 
wave with w Q — * - the beam-confinement factor is 

zero. Even for commonly encountered Gaussian 
beams, this factor is usually very small. For in- 
stance, for \ = 0.5 am and Wq = 50 am. we have s » 
10~ 3 . There is, however, an upper theoretical limit 
to s that is discussed at the end of this subsection. 

The function ib is expanded in powers of s 2 as 

d/ = vii 0 + s 2 \1m + $ 4 \J/ 4 +* . . . . (7) 

The lowest-order term iii 0 represents the fundamental 
mode of the Gaussian beam. By the use of Eq. (5), it 
is easily checked that this mode is 

= iQexp(-iQ(^ 2 + (8) 



Once d/ 0 is known, Eq. to) implies that the higher- 
order functions ^ for n > 1, i.e., corrections to the 
fundamental mode, may be recursively deduced from 




n > 0. (10) 


The functions il/o and are more complicated 1617 
than iJi 0 . Because d/o, . . . depend on ti/ 0 through 


Assume that a Gaussian laser beam is focused to a 
radius wq equal to \ or even \/2, corresponding to s = 
0.16 and $ = 0.32, respectively. The so-called diffrac- 
tion or confinement limit dictates that the beam 
cannot be focused any more tightly than this. The 
existence of this limit may be understood in a qualita- 
tive and intuitive manner as follows. First consider 
a plane wave that is propagating in the positive z 
direction with its electric field polarized in the x 
direction. The variation of E x in the z direction [ i.e., 
one cycle of variation over the distance Az = \i in- 
duces a magnetic field H y . Similarly, the variation of 
H y in the z direction induces a new electric field E x . 
Together the two fields E x and H y recursively induce 
each other, causing the forward propagation of the 
plane wave. Now consider a beam with a Gaussian 
profile in the x-y plane that is again propagating in 
the positive z direction. The additional variation of 
E x in they direction induces a new magnetic field H z . 
and the additional variation of H v in the x direction 
induces a new electric field E z . The variations of E : 
and H : induce yet other fields. Together the new 
fields E z and H, cause the beam to spread transversely 
as it propagates. When E x and H y are slowly varying 
in the x-y plane (i.e., w Q » \ or s 1), the induced 
fields E z and H, are weak and the spreading is slow. 
But when E x and H y are as rapidly varying in the x-y 
plane as they are in the z direction 'i.e., lv 0 « \ or 
5 - 1 2- the induced fields E z and H. are strong, 
and the transverse spreading of the beam is as rapid 
as its forward propagation. The transverse spread- 
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ing resembles the nearly isotropic radiation from a 
point source more than it does a transversely local- 
ized beam that is propagating in a definite and 
unambiguous direction. Therefore the range 
0.16 < s < 0.32 indicates a range for the theoretical 
limit between the directional propagation of a beam 
and an isotropically radiating source. A similar limit 
occurs for Fraunhofer diffraction by an aperture of 
half- width a. For 2a > the diffraction pattern 
contains both relative maxima and relative minima, 
indicating direction dependence. But for 2a <$: \ y 
the diffraction pattern is nearly isotropic in the 
forward hemisphere. 

We may also introduce the £th Davis beam approxi- 
mation defined when only those terms in Eqs. (13)- 
(1S) that explicitly depend on s up to and including 
the power s k are retained. We obtain the first Davis 
beam [k = 1) depending on il> 0 and containing terms 
up to s l ; the third Davis beam (k = 3) depending on 
d> 0 , and containing terms up to s 3 ; the fifth Davis 
beam ( k = 5) depending on iii 0 , i|/ 2 , ii> 4 and containing 
terms up to s°; and so on. 19 None of these beams is 
an exact solution of Maxwell's equations. Maxwell's 
equations are only satisfied in the limit k -+ 


and higher. For k = 3, nonconstant terms occur at 
0(s 8 ) and higher, and they occur at 0(s 12 ) and higher 
for k = 5. Because we understand that the noncon- 
stant terms are artifacts produced by the approxi- 
mate nature of the beam descriptions, they may be 
dismissed. Because the details of the computations 
require much algebra, it is somewhat of a pleasant 
surprise that the resulting g n *s for k = 1, 3, 5 may be 
written in the simple form. 


k 


gn k = 

o 


s 2 * (re - 1)! (n 4- 1 + Q! 
/! (re - 1 - /)! (n + 1)! 


( 21 ) 


which explicitly leads to 

gn • 1 - ( B - lXn + 2Js 2 , (22) 

gn 3 = g« l + l Mn - 2)(n - l)(rt + 2)(n + 3)s 4 
- l /s(n - 3){n - 2)(n - 1 )(n + 2)(re + 3) 
x (« + 4)s 6 , (23) 


gn s = gn 3 + l /24 (re - 4)(n - 3 )(/i - 2 )(/i - l)(/l + 2) 


B. Standard Beams 

In the framework of GLMT, an on-axis laser beam is 
described by the set of BSC’s g n given by 13 

= ‘ \ in ~ l ]jR) ^TTT) f sm2 mf[R ' 9) 

x exp i -iR cos QJP^fcos 0), (19) 

in which r, 0, <b are spherical coordinates (Fig. 1), R = 
kr, j n {R) are spherical Bessel functions, P n l are as- 
sociated Legendre polynomials, and f{R, 0) is defined 
by 

j Er , E 0 \ /cos 6\ 

[if H j = G0S 9)sin e [ sin A * (20) 


x (ft + 3)(rc ■+■ 4 )(n -r 5)s 3 - Vuoin - 5) . . . 
x (n - 1 )(n + 2) . . .(* + 6)s i0 . (24) 

The same procedure may be carried out for the 
general on-axis case with z 0 * 0. The amount of 
algebra is, however, now so great that by-hand compu- 
tations are unreasonable, and use of symbolic compu- 
tation software such as maple is compulsory. Still, 
the g n k ’s end up being given by the relatively simple 
expression 19 


2A+1 

gn“= 2 2 

y- o o 


-2 isz 0 
w Q 


(-lYs* 


(! +j)! 1 

l\j\ l\ 


(re - 1)! (a + 1 + 1)1 

* (n-l-Z)! (» + !)! (25) 


Consider the £th Davis beam approximation and 
denote its radial electric and magnetic fields by E r k 
and H r k , respectively, with k = 1, 3, 5 for the first, 
third, and fifth Davis beams, respectively. These 
lead to the first-order, third-order, and fifth-order 
approximations g n k to the BSC’s in the following 
way. 19 The approximation f k to f is Taylor series 
expanded in powers of the small parameter s, which 
permits an analytical integration of Eq. (19). When 
this is done, it is found that nonconstant terms occur; 
i.e., the result of the integral does not cancel the 
prefactor R/j n (R), contradicting the fact that the 
BSC’s must be constants. The occurrence of such 
nonconstant terms is due to the fact that the approxi- 
mations E r k and H r k do not exactly satisfy the Maxwell 
equations. However, the nonconstant terms appear 
at increasingly higher powers of s when k increases. 
For 6 = 1, the Ois°) and 0[s 2 ) terms are found to be 
constants, with nonconstant terms occurring at Ofs 4 ) 


Equation (25) represents the most general and rigor- 
ous result for the analytical evaluation of the BSC’s, 
which we call the 5 -expansion method. 

Although s is usually small, it is demonstrated in 
Subsection 2.C. that even g rt 3 , which contains terms 
of up to Ofs 10 ), is not sufficient to describe extremely 
focused beams or the BSC’s accurately for large 
partial waves. Because the amount of algebra that 
would be required for higher-order Davis beams to be 
designed and the corresponding g n k> $ to be evaluated 
is extensive, it is appealing to conjecture that Eqs. 
(21) and (25) remain valid for k > 5. The infinite 
generalization then reads 


gn 


;*0 (=0 


= 2 2(-2is^[(-l)'s 2 ' 


UJ 0 


U +j)\ l 
fly! l\ 


in - 1)! in t 1 + /)! 
(re - 1 - /)! in ~ 1)! 


exptikz 0 ). 


(26) 
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which reduces to 


where^,, 1 is given by Eq. (22) and 


A(-l ) l s J !n-l)! (n + l + l)'. „ 

8n ‘ ~ /'• In - 1 - /)! ( n + 1 )! l27 ’ 


forz 0 = 0. 

We call Eqs. (26) and (27) the standard BSC’s, and 
the beam defined by this set of BSC’s is called a 
standard beam. The g n * coefficients of Eqs. (21) and 
(25) will be called the £th-order approximation to the 
standard BSC’s. We claim that standard beams 
should be taken as the ideal description of Gaussian 
beams. This claim will be reinforced by the numeri- 
cal results of Section 3. But before proceeding to 
these numerical results, we must introduce the local- 
ized approximation to the BSC’s of Eq. ( 19). 

C. Localized Approximation 

The localized approximation for a focused Gaussian 
beam 20 is built on the first-order Davis beam of Eqs. 
!8) and (9) and results from the localized interpreta- 
tion of partial- wave expansions. For z 0 = 0, the 
radial electric field of the first-order Davis beam in 
Eq. (11) may be written as 

E r = E q exp( -iAz)sin 0 cos <bf{kr, 0), (28) 


with 


f[kr, 0) = iQ exp - iQ 


r- sin- ( 




(1 - 2Qsr cos 8/ Wq). 

(29) 


The localized approximation g n to the BSC’s g n is 
obtained by application of the localization operator L 
to the function / in Eq. ! 29) according to the prescrip- 
tion 


Lf[R, 0) = fin + y 2 , - 2), (30) 

which is the van de Hulst localization principle 
applied in the focal plane of the beam. 22 The integra- 
tion in Eq. (19) may then be easily performed, 21 
yielding 


g n = exp(-s 2 (n + l / 2 ) 2 ]. (31) 

To motivate the modified localized approximation, 
we now demonstrate that the standard BSC’s of 
Eq. (27) may be approximated by 

g n x = exp[-$ 2 '/r - l)(/i + 2)j. 1 32) 

The demonstration proceeds in the following way. 
The exponential in Eq. (32) may be expanded as 

exp[-s 2 (/i - l)(/i 4- 2)] 

= g n l + Via\ n - 2 ){n - l)|n 4- 2 ){n - 3)s 4 

- V${n - 3)l/z - 2)f n - l>in - 2)(n - 3) 

x in 4- 4 )s 6 + ... f ,33) 


(n - l)(n + 2) 
in - 2 )(n + 3) ’ 


(34) 


0 = 


[n - 1 ) 2 {n + 2) 2 


n - 3)(/i - 2)(rt 4- 3)(n + 4) 


(35) 


The behavior of a and 0 as a function of the partial 
wave n is illustrated in Table 1. For small partial 
waves, when a and 0 are significantly different from 
1, the difference between the exponential in relation 
(32) expanded up to 0(s 6 ) and g n * is small because s 4 
and s s are small (10’ 12 and 10’ 18 for a typical beam 
with s = 10" 3 ). For large partial waves, the 0(s 4 ) 
and 0(s 6 ) terms contribute significantly. But then 
ot = 0 = l with a high accuracy, again validating 
relation (32). The same argument holds for higher 
powers of s as well. 

Relation (32) implies that we may introduce a 
modified localized approximation and a modified local- 
ization operator, 

L mod f[R, 9) = f[(n - l) 1 2 (n + 2) 1 ' 2 , tt/ 21, (36) 
leading to 


gn. mod = exp[-s 2 (n - l)(/i + 2)j. (37) 

The modified localized approximation may also be 
written as 


gn. mod = SXp{-S 2 7irt + Vi) 2 }, 

where 


_ in - l)(>i 4- 2) 

(n + V*f 


(38) 


(39) 


As is shown in Table 1, this ratio also quickly tends to 
1 as n increases. Therefore the localized approxima- 
tion of Eq. (31) is very close to the modified localized 
approximation of Eq. (37). For z 0 * 0, these apprec- 


iable 1 . Coefficients a. $, and 7 of Eqs. i34), 1 35>, and 1 39). 
Respectively, as a Function of Partial Wave 


Partial Wave 

a [Eq. t34)j 

6 i Eq. 1 35 i ( 

7 * Eq. ! 39 ij 

5 

1.166667 

1.314815 

0.925620 

10 

1.038462 

1.144427 

0.979592 

50 

1.001573 

1.005518 

0.999118 

100 

1.0003% 

1.001388 

0.999777 

500 

1.000016 

1.000056 

0.999991 

1000 

1.000004 

1.000014 

0.999998 
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mations generalize to 


-o 


gn = ( 1 4- 2zs — expi ikz 0 ) 
W 0J 


-i 


x exp 


~s 2 {.n ~ V2) 2 


1 *"!" 


^o\’ L 


^«.mod = I 1 + 2is exp(i*z 0 ) 
-s 2 (/i - l)(n + 2) 


x exp 


1 -h 2isz 0/ w 0 


(40} 


(41) 


The results shown in Table 1 indicate that the 
analytical expressions of the localized and the modi- 
fied localized approximations bear a strong resem- 
blance to the infinite-series standard-beam BSC’s. 
But Eqs. (31) and (37) are built on the first-order 
beam, and the standard beams incorporate all the 
higher-order terms. It is a pleasant surprise that 
the localized approximation that is built on the 
first-order beam anticipates these higher-order de- 
scriptions and includes them in an approximate way. 
This will prove to be very useful when the standard- 
beam BSC’s are slowly convergent and is examined in 
more detail in the Section 3. 


3. Numerical Discussion of Localized and 
Standard SSC's 

If we insert the localized and the standard BSC ’3 into 
the beam partial- wave expansions, we generate the 
localized and the standard beams, respectively. We 
could then compare the localized- and the standard- 
beam profiles. This was done for the beam focal 
waist plane in Refs. 19 and 27. In this paper, 
instead, we emphasize the comparison between the 
individual values of the localized and the standard 
BSC’s. 


A. Comparison for z Q = 0 and s 1 

For the case z 0 - 9* we compare the values of the 
BSC’s in Tables 2 and 3 for ( 1) the localized approxi- 
mation of Eq. (31), labeled LA; (2) the modified 
localized approximation of Eq. (37) labeled MLA; and 


(3) the s-ex pansion method for the first-, third-, and 
fifth -order approximations to the standard BSC’s of 
Eqs. (22-24). These are labeled Dl, D3, D5, respec- 
tively. Also compared are (4) the standard-beam 
values obtained from Eq. (21) when k is increased 
until a convergence of 9 significant figures is achieved. 
The values of k for convergence is listed in Tables 2 
and 3, as is the numerical value of g n x . The results 
shown in Table 2 are for the commonly encountered 
situation of s = 0.001. 

First we consider the convergence of Eq. (21) for 
the standard-beam BSC’s for s = 0.001, which are 
typical of laser Doppler and phase Doppler instru- 
ments, 1 ” 3 and that correspond to focusing of the beam 
to a radius u/ 0 - 150\. The evaluation of g n k has 
been carried out by the use of the symbolic computa- 
tion software maple. This is compulsory because 
when one is evaluating the standard coefficients, the 
number of digits required in the computations to 
obtain 9 significant figures in the results may be far 
beyond what is available with fortran double preci- 
sion variables. MAPLE allows one to carry out evalua- 
tions with an arbitrary number of significant figures, 
which is only limited by the host-computer available 
storage, by setting the maple variable digits to a 
prescribed value. For instance, digits = 12 is enough 
to evaluate g n k for small n. For n = 2500, the 
evaluation of g n 31 requires digits = 20, and for n = 
5000, the evaluation of g„ 101 requires digits = 40. 
Therefore, although the standard BSC’s provide 
benchmark values for the BSC’s that describe a 
focused Gaussian beam, such benchmark values may 
in practice be difficult to obtain for large partial waves 
and tightly confined beams. 

Examining the sequence ofg^’s for various partial 
waves n and various beam orders k y we may follow the 
convergence of the standard scheme. Up to the 
partial wave n = 5, g n l is sufficient; i.e., the standard 
BSC’s are correctly evaluated only by the use of 
first-order Davis beam. For 10 < a < 100, the use 
of the third-order Davis beam is required. Eventu- 
ally it is necessary to rely on the Ath-order standard 
scheme with k values larger than 5. For example, 
for n = 1000, 2500, and 5000 we need k = 15, 31, and 
101, respectively. Large partial waves n are assoei- 


Table 2. BSC’s as a Function of Partial Wavs for $ = 0.001 for the Localized Approximation i LA); the Modified Localized Approximation i MLA); the 
First- * 01 ), Third- i 031, and Fifth-order 1 05) approximations to the Standard Beam; and the Standard Beam < S)‘ 


n 

LA 

MLA 

Dl 

D3 

D5 

*,S 

1 

0.999997750 

1.000000000 

1.000000000 

1.000000000 

1.000000000 

l, Same as Dl 

2 

0.999993750 

0.999996000 

0.999996000 

0.999996000 

0.999996000 

1, Same as Dl 

5 

0.999969750 

0.999972000 

0.999972000 

0.999972000 

0.999972000 

1, Same as Dl 

10 

0.999889756 

0.999892006 

0.999892000 

0.999892006 

0.999892006 

3, Same as D3 

50 

0.997452999 

0.997455243 

0.999452000 

0.999455238 

0.997455238 

3, Same as D3 

100 

0.989950586 

0.989952813 

0.080902000 

0.989952793 

0.989952793 

3. Same as D3 

1000 

0.367511653 

0.367512480 

<0 

0.332834669 

0.366292083 

15, 0.367511867 

2500 

0.001925633 

0.001925638 

<0 

<0 

<0 

31, 0.001925639 

5000 

0.138186 x lO" 10 

0.138X87 x lO’ 10 

<0 

<0 

<0 

101. 0.138208 x lO' 10 


°For the standard beam, the number of terms in the infinite senes of Eq. <27} required for convergence to 9 significant figures \k\ is also 
given. 
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Table 3. 

SSC's as a Function of Partial Wava for a - 0.16 for tha Localized Approximation (LA); the Modified localized Approximation ( MUti- the 
First- 1 01 ), Third- 1 D3), and Fifth-Order I D5) Approximations to the Standard Beam; and the Standard Beam (S1« 

n 

LA 

MLA 

Dl 

D3 

D5 

k,S 

1 

2 

4 

6 

10 

15 

20 

25 

0.944027482 
0.852143789 
0.595472542 
0.339052607 
0.059463060 
0.002132629 
0.000021266 
0.589603 x L0- 7 

1.000000000 
0.902668412 
0.630778820 
0.359155441 
0.062988600 
0.002259075 
0.000022526 
0.624562 x 10" 7 

1.000000000 

0.897600000 

0.539200000 

<0 

<0 

<0 

<0 

<0 

1.000000000 

0.897600000 

0.616138215 

0.327063245 

<0 

<0 

<0 

<0 

1.000000000 

0.897600000 

0.618138215 

0.343026339 

<0 

<0 

<0 

<0 

1. Same as Ol 
1, Same as Ol 
3. Same as D3 
5, Same as D5 
9, 0.058365667 
15, 0.002267813 
19, 0.000031912 
25, 1.853835 x 10' 7 


'For the standard beam, the number of terms in the infinite senes of Eq. (27) required for convergence to 9 significant figures Ik) is also 


ated through the localized interpretation with geomet- 
ric light rays that are passing far from the beam 
axis. 22 * 23 Therefore the description of the outer parts 
of the beam requires higher k orders than the descrip- 
tion of the central region. This observation must be 
reconciled with the previously demonstrated fact 18 
that a first-order Davis beam satisfies MaxwelTs 
equations up to Ofs 2 ) uniformly over all space. Our 
results on the convergence of the standard scheme in 
Table 2 indicate that the situation is more subtle 
because the coefficients a* of the various powers of s 
can make terms such as k > 2, significant if a* is 
an increasing function of the partial wave n and if n is 
big enough. Clearly, for the on-axis case, geometric 
rays associated with large partial waves possess van- 
ishingly small amplitudes that are ineffective in the 
light-scattering process, so that a poor evaluation of 
the corresponding BSC’s should not be influential. 
Note, however, that whether a partial wave is effec- 
tive also depends on the size of the target particle 
through the Lorenz-Mie partial-wave scattering am- 
plitudes a n and b n . In addition, low partial waves are 
classically associated with backscattenng, and large 
partial waves with side scattering. Thus, when one 
compares g n values such as those in Table 2, where 
gi = 1.0andg 500 o 3:5 10" n , which is vanishingly small 
in comparison with g u we should actually compare 
the light scattered in different directions. A more 
refined discussion should then take into consider- 
ation scattering diagrams in an actual scattering 
process. 

We now consider the accuracy of the localized 
approximation and the modified localization approxi- 
mation BSC’s when compared against the benchmark 
standard-beam BSC’s for s = 0.001. In Table 2. the 
comparison between the modified localized approxima- 
tion and the standard scheme is excellent. Up to the 
partial wave n = 100, the difference between the 
modified localized approximation and the standard 
scheme typically does not exceed 1 part in 10 8 . Even 
for n = 5000, the disagreement lies in the fifth 
significant figure. There the modified localized ap- 
proximation based on the first-order Davis beam 
anticipates the information contained in the 101st 
order of the standard scheme. This unexpected inter- 
nal coherence is considered as a cross-check of the 


validity of the modified localized approximation, and 
of the fact that standard beams should be considered 
as the ideal reference beams. Finally, the localized 
approximation agrees reasonably well with the stan- 
dard BSC values. But the agreement for the modi- 
fied localized approximation is better, especially for 
n < 100. 

8. Comparison for z 0 = 0 and s =0.16 

Table 3 now provides a comparison for z 0 = 0 and s = 
0.16 [i.e., 1/(2—)] near the theoretical confinement 
limit. The range of important partial waves is much 
smaller than in Table 2. This is a direct consequence 
of the localized interpretation; i.e., a BSC of partial 
wave n is associated with the amplitude of the 
geometric light ray that is passing at a distance 

{n + u>)\ 

P* = ~ 2 - (42) 

from the beam axis at the focal waist. From relation 
(32) it can be seen that the amplitude decreases to 
1/e 2 of its value on the z axis for n = 1500 if s = 10“ 3 
and for n ~ 10 if s = 0.16. The n values are strongly 
correlated with the necessity of using bigger k orders 
to obtain convergence for the standard-beam BSC’s. 

There is also an increase in the difference between 
the modified localized approximation and the stan- 
dard values. This difference is dramatic for n — 25. 
There is also an increased difference between the 
localized and the modified localized approximations. 
Also, depending on the partial wave n , the modified 
localized approximation may compare more favorably 
or less favorably with the standard BSC values than 
the localized approximation. The deterioration of 
the comparisons for s = 0.16 is consistent with the 
approach to the physical confinement limit. 

C. Comparison for Zq * 0 

We now focus our attention on the convergence of the 
standard-beam BSC’s for the general on-axis case 
z 0 ~ 0, withg** given by Eq. (26). Computations are 
carried out by means of a maple procedure, increasing 
k in Eq. 1 25 j until the convergence test g n * = g n k ~ 2 is 
satisfied to an accuracy of 50 significant figures. 
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The value of k for which convergence is reached is 
denoted by K . In Fig. 2, K is displayed versus z Q for 
gi x , the BSC for the first partial wave. The beam- 
waist radius u/ 0 is used as a parameter. The number 
of terms required for convergence increases when z Q 
increases and when w Q decreases, i.e., when the focal 
waist of a tightly focused beam is far upstream or 
downstream from the origin of coordinates. The 
increase versus z 0 is particularly sharp for the most 
focused beam (u/ 0 = 0.25 jim, s = 0.32) at the upper 
limit of the physical confinement range. Figure 3 
presents the same data shown in Fig. 2 versus the 
dimensionless quantity z Q /L, in which l is the spread- 
ing length. Because of the fact that l is the natural 
characteristic length to rescale the 2 coordinate, all 
the curves in Fig. 2 collapse to a single curve in Fig. 3. 

In Figs. 2 and 3 only the first partial wave was 
considered. To extend the analysis to all partial 
waves, K is presented in Fig. 4 as K (D t z 0 /Z). In Fig. 
4 z 0 is still rescaled by l. Rather than the partial 
wave n, the ordinate is now taken to be 

D = (n - l) 1/2 (n + 2) 1/2 s, (43) 

which is the distance p from the beam axis associated 
with n through the modified localized interpretation 
[Eq. (44), below] and rescaled by u/ 0 . Figure 4, for 
X = 0.5 |xm, wq = 5 p.m, and s - 0.016, demonstrates 
how the number of terms required for convergence 
increases when z 0 or n increases. 

To some extent, these results might be considered 
troubling. Consider, for instance, a small value for 
s, which leads us to expect that the beam may be 
safely described by a first-order Davis beam. Figure 
4, however, tells us that this conclusion is true only in 
a small region that surrounds the beam-waist center 
and that the standard-beam BSC’s will be slowly 
convergent otherwise. Fortunately, as mentioned 
above, the localized approximation anticipates the 
behavior of high-order Davis beams and therefore 
may be used in this case as an alternative to the 


order k versus Zo, for g, 



Fig. 2. Value K of the £th-order standard beam required for 
convergence of the BSC g\* as a function of z$. 



Zo/l 


Fig. 3. Value K of the fcth-order standard beam required for the 
convergence of the BSC gi" as a function of zq/1. The individual 
curves from Fig. 2 now coincide. 


standard-beam procedure to provide a fast and accu- 
rate way to evaluate BSC’s. This is the most impor- 
tant result of this paper. When one is computing 
light scattering with the GLMT formalism, the de- 
scription of the incident beam should be accurate and 
should permit scattering calculations to be performed 
rapidly. The standard-beam BSC’s yield the best 
description of the beam. But in certain circum- 
stances their slow convergence causes a large increase 
in the computer run time of scattering calculations. 
The localized approximation, on the other hand, 
satisfies both criteria of accuracy and computational 
speed, making it a useful tool in GLMT calculations. 

The comparison between the standard, localized, 
and modified localized BSC’s may be complicated in 


0.00 0.10 0.20 0.30 0.40 



z«/l 

Fig. 4. Value K of the fcth-order standard beam required for the 
convergence of the BSC g n x as a function of z 0/ 1 and D = 
sin - l) 1 + 2) 1 2 for X = 0.5 u.m. ll/q - 5 y.m. and s - 0.016. 
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the so-called off-axis case, when only one side of the 
particle is illuminated by the beam. To investigate 
this problem, it is necessary to design an expansion 
method for higher-order Davis beams and correspond- 
ing off-axis standard-beam expressions. As of yet, 
this task has not been undertaken, to our knowledge. 

4. Top-Hat Beams 

Originally, the localized approximation for a focused 
Gaussian beam was obtained by analogy to van de 
Hulst’s localization principle for plane waves. This 
led us to associate a BSC g n with the amplitude of a 
geometric light ray that is passing at a distance p n 
from the beam axis, in which p* is given by Eq. (42). 
In the modified localized approximation, we have 

Pn - (n - 1 ) l 2 (n + 2) 1 2 ~ (44) 

These relations allow us to investigate light scatter- 
ing by the use of the localized approximation for 
beams more general than Gaussian beams. For 
instance, consider top-hat beams, which are used in 
certain optical particle-sizing instruments. 7 *" 10 Top- 
hat beam scattering has been calculated previously 
with the localized approximation, without, however, 
directly demonstration of the validity of the localized 
approach to the case of top-hat beams. 10 An assess- 
ment of this validity is provided here, illustrating the 
potentialities of the localized interpretation. We con- 
sider an idealized beam profile defined by Eq. (20) in 
which 


H in the x-y plane. The explicit forms of the func- 
tions F u F 2 , F 3 in terms of the BSC’s are given in 
Ref. 19. 

For the localized beam model of Eq. (46) the 
dominant shape function F x [kr) was calculated for the 
top-hat profile laser beam, for w Q = 25 p.m, 7.5 pm, 
and 2.5 pun, and is shown in Figs. 5(a)— 5(c). The 
results are generally encouraging. The fields are 
virtually constant from the z axis out to the radius w Q 
in the x-y plane as hoped. But beyond w 0t instead of 
being rigorously zero, the fields are oscillatory, with a 
slowly decreasing amplitude. This is reminiscent of 
the oscillatory ringing in the Fourier transform of a 
function with a hard edge. 29 As the radius of the 
top-hat beam w 0 decreases, the oscillations become 
coarser and their amplitude increases. F or example, 
when w 0 = 25 pm, the amplitude of the oscillations 
falls by an order of magnitude from its value in the 
plateau region when the distance p from the z axis is 
approximately 1.2 w 0 . For w 0 = 2.5 pm this occurs 
when p =» 1.7 w 0 . The exact effect that these oscilla- 
tions in the localized electric and magnetic fields have 
on the far-field scattered intensity calculated in the 
GLMT framework is not known. But it is expected 
to be small because the amplitude of the oscillations is 
only a small fraction of E 0 , the field strength in the 
plateau region. 

The sharp cutoff in Eq. (45), however, is not 
observed in oscilloscope traces of experimental beam 
profiles. 7 In actuality, the fields possess a smooth 
but rapid roll-off for p > w 0 . A more realistic model 
of a top-hat beam is then 


fikr, 9) = 


if r sin 9 < u/ 0 
if r sin 9 > u/ 0 , 


(45) 


fikr, 9) = 


1 

exp[-(r sin 9 - w 0 ) 2 /*r] 


if r sin 9 ^ uj q 
if r sin 9 > w 0 , 


corresponding to the plateau region of constant illumi- 
nation of radius w 0 around the beam axis. If we rely 
on the localized interpretation of Eq. (30), such a 
beam should be generated by localized BSC’s given by 


(50) 

where e is the small roll-off distance of the beam in 
the x-y plane. The localized beam model for Eq. { 50 ) 
with the prescription of Eq. (30) is then 


gn 


1 if n <> kw 0 - l /2 
0 if n > kw 0 - l h 


(46) 


The accuracy to which Eq. (46) describes a top-hat 
beam was tested in the following way. It is known 
that a given set of BSC’s g n defines an on-axis laser 
beam that is an exact solution of Maxwell’s equations. 
The electnc-field components of this beam in the 
beam- waist plane (i.e., 9 = - 21 are given by 

E x {kr , tt/2, 6/ = F x {kr * - F 2 {kr)sin z 6, (47) 

E y {kr, 77/2, d>) = F 2 (kr\sin d> cos <b, (48) 

E z ikr, 77 2, (b) = F 3 {kr) cos d>, (49) 

with similar relations for the magnetic fields’ compo- 
nents. The function F { describes the dominant shape 
of the beam profile, and F- z and F z denote corrections 
to the dominant shape induced by variations in E and 


gn “* 


exp(-(n + V 2 - kwqf/kH 2 } 


if /l ^ kW(y - l /o 
ifn > kw 0 - V 2 . 


(51) 


The dominant shape function F x of the beam defined 
by Eq. (51) was calculated for w 0 = 25 p.m and e = 
0.05ii/ 0 , and is shown in Fig. 6(a). The oscillations in 
the field persist for p > w 0 . But their amplitude has 
now been decreased to less than 10“ * E Q , and their 
effect on the far-field scattered intensity is similarly 
reduced. This decrease is confirmed in Fig. 6(b) for 
wo = 25 p.m with the more gradual roll-off € = 0.1u/ 0 . 
The oscillation level has now decreased to less than 
10“ 4 E o amd is again reminiscent of the reduction in 
the oscillatory ringing in the Fourier transform of a 
function with a smooth, gradual edge. 29 The results 
in this section therefore further support the localized 
interpretation of partial-wave analyses and provide a 
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Fig. 5. Dominant beam-shape function Fi for the localized top- 
hat beam of Eq. (45) with iai ujq =* 25 jxm, (b) wq * 7.5 p.m, and (c) 
* 2.5 ^tm as a function of the distance p from the z axis. 




ibi 


Fig. 6. Dominant beam-shape function F\ for the localized top- 
hat beam of Eq. ( 50 ) T which possesses a smooth roll-off of width 
e. The curves are for o/ 0 = 25 p.m and la) « = 0.05u;o and (bl e - 
O.lOi^o. The smoothing of the edge of the beam profile dramati- 
cally decreases the oscillatory ringing in F\. 


new localized description of top-hat beams [eq. (51)], 
which improves on the description of Eq. (46) previ- 
ously published in Ref. 10. 

5. Conciusion 

In the GLMT framework, the question what are the 
electromagnetic fields of a laser beam is equivalent to 
the question what are the values of the BSC's that 
describe the beam? By the use of a so-called s-ex- 
pansion method, one may obtain the fields of stan- 
dard beams associated with standard BSC’s. We 
claim that these coefficients represent an ideal descrip- 
tion of a Gaussian beam. Standard beams make 
possible the study of both mildly focused and ex- 
tremely focused beams. The so-called localized ap- 
proximation received a rigorous justification. It pro- 
vides a simple analytical expression for the BSC’s 
that is very close to standard beam values. Finally, 
the localized interpretation of partial-wave expan- 
sions that underlies the localized approximation re- 
ceived further support from the investigation of 
top-hat beams, because they may be accurately de- 
scribed with the localized approximation. 

This work was supported in part by National 
Aeronautics and Space Administration grant NCC-3- 
204. 
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